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The compressible, turbulent, time dependent and three dimensional flow in a 
labyrinth seal can be described by the Navier-Stokes equations in conjunction with 
a turbulence model. Additionally, equations for mass and energy conservation and an 
equation of state are required. To solve these equations, a perturbation analysis 
is performed yielding zeroth order equations for centric shaft position and first 
order equations describing the flow field for small motions around the seal center. 
For numerical solution a finite difference method is applied to the zeroth and 
first order equations resulting in leakage and dynamic seal coefficients respec- 
tively. 
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Forces on the shaft in z, y direction 

shaft displacements 

direct and cross-coupled stiffness 

direct and cross-coupled damping 

axial, radial and circumferential velocity 

pressure, turbulence energy, turbulent energy dissipation 

temperature 

laminar, turbulent and effective viscosity 

(p e = 

density 

specific heat for constant pressure 
time 

axial, radial and circumferential coordinate 
radial coordinate after transformation 
turbulent Prandtl number 
constants of the k-e model 

constants of the k-e model 

production term in turbulence model 
laminar dissipation 
general variable 
general source term 

seal clearance for centric shaft position 

seal clearance for eccentric shaft position 
radius of the precession motion of the shaft 
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e = r Q /C o perturbation parameter 

o rotational frequency of the shaft 

fl precession frequency of the shaft 

Tj radius from seal center to end of seal fin 

r radius from seal center to stator 

a 

Transformation constants 


SUBSCRIPTS 

0 zeroth order variables 

1 first order variables 

R Rotor 

S Stator 

c cosine 

s sine 


INTRODUCTION 

The problem of subsynchronous vibrations in high performance turbomachinery has 
been investigated for many years. However, for some destabilizing effects there is 
no satisfactory solution from the qualitative as well as quantitative points of 
view. In machines like turbocompressors and turbines, gas seals and especially 
labyrinths have an Important influence on the dynamic behavior. These elements with 
turbulent flow conditions have the potential to develop significant forces 



which may lead to self-excited vibrations of the shaft. 

To predict the stability behavior of a high performance rotating machinery, the 
rotordynamic coefficients (eq. 1) of the turbulent seals have to be known. Many 
theoretical studies, numerical calculations and measurements have been carried out 
to determine the stiffness and damping characteristics. Most of the recent theore- 
tical developments are based on a bulk flow theory, considering the shear stresses 
only at the wall but not within the fluid. Nelson (ref. 1) has derived a procedure 
with a ’one volume* model for straight and tapered seals. Iwatsubo (ref. 2) as well 
as Childs (ref. 3) have also used the one volume model for labyrinth seals. How- 
ever, in case of such complicated seal geometries difficulties may occur, when 
neglecting the local turbulence characteristics within the fluid, particularly at 
locations with strong velocity changes. 

Wyssmann et al. and Wyssmann (ref. 4 and 5) have extended the bulk flow theory for 
labyrinth seals, working with a two volume model. The latter takes into account the 
shear stresses between the core flow and the Jet flow. Wyssmann* s two volume model 
is derived from results of a finite difference calculation for a rotationally sym- 
metric single cavity turbulent flow (centered position of the shaft), based on the 
time averaged Navier-Stokes equations with a k-e model. A first attempt to apply 
finite difference methods to labyrinth seals was made by Rhode et.al. (ref. 6), who 
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calculated only one chamber seals for centric shaft position. However, bulk flow 
theories show a substantial lack In the prediction of labyrinth seal coefficients. 
As already shown by Nordmann, Dletzen and Weiser (ref. 7), the solution of the 
time-averaged conservation equations for momentum, mass and energy in conjunction 
with a k-e turbulence model by means of a finite difference technique is a physi- 
cally more realistic way for calculating rotordynamlc coefficients of gas seals. 

This paper Is concerned with the extension of the already mentioned calculation 
procedure for annular seals to straight-through labyrinths. 


GOVERNING EQUATIONS 

To describe the flow in a labyrinth seal, we use the time-averaged conservation 
equations for momentum, mass and energy and the equation of state for a perfect 
gas. The correlation terms of the turbulent fluctuation quantities are modelled via 
the k-e turbulence model of Launder and Spalding (ref. 8). 

All these differential equations can be arranged in the following generalized form: 

I E< p *) + r §F< prv *) + r §e< pw *) + §£< pu *) = 

( 2 ) 

I 3_/ rr + I «_rr 1 s 

r dr* dr } r d0* J r dQ' dx*' dx' <f» 

where ^ stands for any of the dependent variables, and the corresponding values 
of r<p and S+ are indicated in Table 1 (see Appendix A). 


PERTURBATION ANALYSIS 

Obviously, this set of equations is not solvable due to the limited performance 
capabilities of today's computers. Therefore, a stepwise procedure is performed in 
order to reduce the three dimensional and time dependent flow problem to a two 
dimensional time independent one. 

First, the governing equations are transformed to another coordinate system by 
introducing a new radial coordinate q whereby the eccentrically moving shaft is 
converted to a shaft rotating in the center of the seal (Fig. 1). 


Fig. l: 
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The transform functions are given In Fig. 2. 



Strips on rotor: 

C 6(r-r, ) 
o v 1' 


rj=r + 


1 (C^Hr-rjJ+Cfi 



Fig. 2: Transform functions 


Assuming small shaft motions on a circular orbit around the centric position we 
perform a perturbation analysis by Introducing the following expressions for the 
dependent variables into the governing equations: 

u = u+ eu. v = v+ ev, w = w + e w, 
o 1 o 1 o 1 

p = P o + e pj T = T q + e Tj p = p q + e p { (4) 

This procedure yields a set of zeroth order equations governing the centric flow 
field and a set of first order equations describing the flow field for small 
eccentric shaft motions. 

The first order equations depend also on the circumferential coordinate 0 and on 
time t. To eliminate the circumferential derivations, we insert circumferentially 
periodic solutions for the first order variables in the corresponding set of equa- 
tions: 
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After separating the equations into sine and cosine terms and rearranging them by 
introducing the following complex variables. 
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( 6 ) 


U 1 = u lc * lu ls V 1 ' v lc * lv ls *1 c *lc * lw U 

f l * p lc * lp ls T 1 “ T lc + 1T ls p l ‘ "lo * lp l s 

a circular shaft precession orbit (Fig. 3) with frequency fi and corresponding 
solutions for the first order variables are assumed whereby the temporal deri- 
vations can be eliminated: 



Finally the zeroth and first order equations are arranged to a generalized form'- 


d .. d, 9*. 1 a, . 1 d, r 

9x^ P o U o^ $ dx ^ T) d^oV* q 3q^ <f> ^ dr^ 


* 


( 8 ) 


Table 2 (see Appendix B) indicates the values for the zeroth and first order 
equations. 

The constants Dj-D^q (given in the Appendix C) arise from the coordinate trans- 
formation and depend only on zeroth order variables. These terms are zero for the 
labyrinth chambers. 


FINITE DIFFERENCE METHOD 

For the numerical solution of the governing equations we apply a finite difference 
technique in the same manner as suggested in many publications (see e.g. ref. 9). 
First the calculation domain is covered by a finite difference grid (Fig. 4). Dis- 
cretizing the equations results in algebraic expressions linking every node point 
to his four neighbours. Solution of zeroth order equations results in flow field, 
pressure distribution, and leakage. For first order solution, the equations for the 
turbulence quantities k, e can be dropped due to assumed small rotor motions. 
Because Eq. (1) contains four unknown parameters (K.k.C.c). first order solution is 
performed for two different shaft precession frequencies 0. This results in two 
sets of dynamic forces Fj , F^ obtained from the integration of the pressure pertur- 
bations (see Eq. 9). 
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( 9 ) 



Finally the dynamic coefficients are calculated from Eq. 10. 

F , 

- ~ = K + cO 

o 

F 2 

- — = -k + cn (io) 
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Fig. 4: Nonuniform finite difference grid for labyrinths 


BOUNDARY CONDITIONS 


- Zeroth order equations 

For the velocities the near wall regions are represented by the logarithmic wall 
law. The entrance conditions are iteratively calculated depending on entrance 
Mach number and loss coefficient. The entrance swirl is assumed to be known. 

- First order equations 

Due to shaft precession and coordinate transformation, one has to pay attention 
to the following boundary conditions for the radial and circumferential velo- 


cities: 


v jg = (0* • 0. ) 

A 

W 1S = (°* * 0 *) 

V 1R = (°* * ( n-w ) c 0 ) 

W 1R = ( n V°-> 
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For the pressure, the exit perturbation is zero in case of exit Mach number less 
than one ("undercri tical” flow) which occurred in the appl i cat ion examples. 


APPLICATION 

For testing our program, we made comparative calculations to measurements performed 
by Benckert (ref. 10). 


Seal Data: 

Shaf t diameter 
Strip height 
Strip pitch 
Clearance 


0, 150 m 
5,5 mm 
8 mm 
0,5 mm 


Number of chambers : 
Entrance pressure p^ : 

Exit pressure p & : 

Reservoir temperature: 

Laminar viscosity 
Perfect gas constant : 
Specific heat ratio : 
Rotor surface velocity: 


3 

1,-425 bar 
0,950 bar 
300 K 

1.8-10" 5 Pa s 
287.06 J/kgK 
1.4 
0 m/s 


The calculation grid, the flow field and the pressure distribution are shown in 
Fig. 5-7. 


In Fig. 8-10 mass flow, direct and crosscoupled stiffness are plotted versus the 

9 

swirl parameter, as defined by Benckert (E = 0,5 *p *w / (p -p )). 

v n r n a vr n r a * J 



Fig. 5: Finite difference grid 
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Fig. 8: Hass flow 
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Fig. 9: Direct Stiffness 
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Fig. 10: Crosscoupled Stiffness El [ ] 


CONCLUSION 

We have shown that a finite difference method based on the Navier-Stokes equations 
in conjunction with a turbulence model is more suitable to describe the dynamic 
coefficients of gas seals than the simpler models based on bulk flow theories. Of 
course the calculation time is also much greater. But great progress is still 
possible. For example, advanced solution algorithms like Multi-Grid, Newton-Sche- 
mes, SIP (strongly implicit procedure) can be implemented to speed up the computer 
program. 

Further on. a three dimensional calculation procedure is currently under develop- 
ment to determine the coefficients of more complex geometries like stepped and 
interlocking labyrinths (Fig. 11,12). 



Fig. 11: Stepped seal 
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Production term in the equations for the turbulence quantities: 

c = u r2f f— i 2 +f— ) 2 +/I a V^ 2 m aY t a< h 2 ir 1 av t aw _ ii 2 

M e L - 4 ll ar^ *'-dx ) n r aSr 1 '^a^aF' +l r a&dF P 
,dw l 3u.2 2 ,du 1 dw l 8 , ,,2, 

+( aF + 7 30^ ' 3 <aF + Fae + FaF < rv » 3 

Laminar Dissipation in the energy equation: 
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APPENDIX B 

Table 2a: Zeroth order equations 



Equation of State: 

o 


Table 2b: First order equations 


Equation of State: 
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Only the first order continuity-equation to determine shows a slightly modified 
form. 
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APPENDIX C 


First order transformation constants 

Note-* r = r for strips on stator 
x a 

r = r. for strips on rotor 
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